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Abstract 

Context. Direct numerical simulations of plasmas have shown that the dynamo effect is efficient even at low Prandtl numbers, i.e., 
the critical magnetic Reynolds number Rm f that is necessary for a dynamo to be efficient becomes smaller than the hydrodynamic 
Reynolds number Re when Re — » oo. 

Aims. We test the conjecture that Rm c tends to a finite value when Re — » oo, and we study the behavior of the dynamo growth factor 
y at very low and high magnetic Prandtl numbers. 

Methods. We use local and nonlocal shell models of magnetohydrodynamic (MHD) turbulence with parameters covering a much 
wider range of Reynolds numbers than direct numerical simulations, that is of astrophysical relevance. 

Results. We confirm that Rm c tends to a finite value when Re — > oo. As Rm — > oo, the limit to the dynamo growth factor y in the 
kinematic regime follows Re^, and, similarly, the limit for Re — > oo of y behaves like Rrrr 9 , with x p' « 0.4. 
Conclusions. Our comparison with a phenomenology based on an intermittent small-scale turbulent dynamo, together with the 
differences between the growth rates in the different local and nonlocal models, indicate that nonlocal terms contribute weakly to the 
dynamo effect. 
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1. Introduction 

The magnetic Prandtl number Pm (the ratio of the kinematic 
viscosity v to the magnetic diffusivity rj) is one of the non- 
dimensional parameters that control the properties of plasma. 
This ratio translates to Pm = Rm/Re, where Re and Rm are the 
hydrodynamic and the magnetic Reynolds numbers respectively, 
and this means in particular that it is an important parameter in- 
fluencing the properties of turbulence in a plasma. Its value is 
large in hot, rarefied plasmas such as the interstellar or intraclus- 
ter medium and the solar wind, and it is small in cool and dense 
plasmas as in planets and the Sun's convective zone, which is at 
the origin of the solar magnetic field. 

|IskakovetaTlP007| l and |Sche kochihi n et al.] ( |2007| > demon- 
strated with direct numerical simulations that the turbulent dy- 
namo is effective even for small Pm, but the minimum Prandtl 
number that they were able to use was limited by the resolution 
of the simulation to 0.070 (with eighth-order hyperviscosity). 
Owing to this limitation, they cannot tell for sure whether the 
critical magnetic Reynolds number Rm e tends to a finite value 
when the hydrodynamic Reynolds number Re tends to infin- 
ity (their maximum Re is 6200). Large- scale dynamos at Pm 
down to 10" 3 have also been obtained by Brandenburg (2009 1. 
However, the range of Pm that is relevant to astrophysics is even 
wider. 

shell models of magnetohydrodynamical (MHD) turbulence 
allow us to go much beyond these limits, up to Re « 10 12 
and Pm ss 10 ±I2 in this paper. They are dynamical models 
of the nonlinear interactions between fields on different scales 
and have been developed for the study of turbulence in vari- 
ous frameworks such as hydrodynamics ( Gledzer|1973 1, MHD 



(Nigr o et al.|2004l [Buchlin & Velli|2007| l, Hall-MHD (Ga ltier& 
Buchlin||2007 1, and MHD with a global rotation rate (|Perrone 
et al.|201 l| l. A review of shell models can be found in |Biferale 



( |2003| > and in the book by |Bohr et aE1 ( |2005] l. 

The simplification that they provide, making them comple- 
mentary to direct numerical simulations (that are limited by the 
resolution one can afford), has allowed numerous results to be 
obtained. MHD shell models have been proven to display dy- 
namo action (|Gloaguen et al.|[l985| |Frick & Sokoloff|[ l998; 
Sahoo et al. 2010), including magnetic field rever sals (|Perrone 
et al.||201 \\. T hey have in particular been used by Stepanov & 



|Pluni an (2006) for the dynamo at low Pm, and later by the same 
authors ( Stepanov & Plunian 2008 ) to propose a phenomenology 
of the turbulent dynamo at both low and high Pm (using nonlocal 
shell models: Plu nian & Stepanov|2007| l. 

In this paper, we use local and nonlocal shell models to 
derive properties of the dynamo over a very wide range of 



Reynolds and Prandtl numbers, extending the results of Iskakov 
etal](|2007| and |Stepanov & Plunian| ( |20T)6] i. 



UGloaguen et al.|1985[|Yamada & Ohkitani|1987| l, reduced MHD 



2. Model equations and numerical set up 

In shell models, the Fourier space for the fields of MHD is di- 
vided in concentric shells of radii k n - A" (A is the separation 
factor between shells), and the fields in each of these shells 
are represented by the complex scalars u„ for the velocity field 
and b„ = B„j yj^iop for the magnetic field (p is the density of 
the plasma, b„ has the dimension of a velocity). The nonlinear 
terms of incompressible MHD, a convolution in Fourier space, 
are written in the following symmetric form when expressed as 
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a function of the Elsasser variables z* = u n ± b n 

( d ^t) NL = l ^-{Q n {z\z\a-b) + Q n {z*,z\a + b)). (1) 

In general, the coefficients of the nonlinear terms are determined 
by the conservation of quantities representing the invariants of 
incompressible MHD in three dimensions (3D), namely the en- 
ergy, cross helicity, and magnetic helicity. 

In the local "GOY" she ll model ( |Giuliani & Carbone|1998) 
Stepanov & Plunian 2006), the nonlinear terms are limited to 



quasi-local interactions between three consecutive shells 

Onix, y c) = Cl x; +l y; 1+2 + c 2 r n _ x r n+ , + c 3 r n _ 2 r n _, (2) 

with c - a + b and 

fll = 1 a 2 =(l-A)/A 2 a 3 = -l/A 2 
bi=b 2 =b 3 = l/A(A+ 1). 



In the nonlocal "Sabra" model ( Plunian & Stepanov]|2007 
Stepanov & Plunian 2008 ), where the hypothesis of quasi-loca 
interactions is released, 



N 

Qn(X, Y, c) = 2^ T m (cj n X n+m Y„ +m+ i + c 2 n X n 

m=\ 



~*~ c m-Xn-m-l Y n -l) 



(3) 



with the coefficients 



a l m = A m (A+\) ai = -A-(-AT 



= (1 - (-AT m )l A 



1 = 1 

m 1 

JA(A- 



1) 



and where N is the number of shells used in the computation. 
The a exponent in T m controls the strength of long-range (in 



Fourier space) nonlinear interactions: following |Stepanov & 
Plunian (2008), we use a — -1 for strong nonlocal interactions, 
and or = -5/2 for weak nonlocal interactions; a - -oo would 
correspond to no nonlocal interactions, i.e., to a local model. 

The forcing /„ is a solution of a stochastic Langevin equation 
and is applied to the two first modes n — and 1 of the veloc- 
ity only, and the dissipation is modeled by Laplacian diffusivity 
coefficients v (viscosity) and rj (magnetic diffusivity). 

The equations for all 3 models (local GOY model, weakly 
and strongly nonlocal Sabra models) are solved numerically us- 
ing a 3rd-order Runge-Kutta scheme for the nonlinear and forc- 
ing terms, and a first-order implicit scheme for the dissipa- 
tion. The time step is adaptive and is set to the smallest time 
scale of the nonlinear terms, with a security factor of 5. We 
use A — 2, and the initial condition for the velocity field is 
u„ = iMole'^^V*' , where \u \ = 10 _1 , St = 10~ 3 , and <f>„ 
are independent random phases. After one large eddy turnover 
time £j with no magnetic field (the shell-model equations reduce 
to hydrodynamic shell-model equations), we introduce an ini- 
tial magnetic field b„ = |6 |e ! '< ,> »^ 1/ V''*>°, where |i | = 10" 10 
and (p' n are also independent random phases: the cross helicity is 
close to zero, and the Lorentz force (scaling as k 2 ub) is negligi- 
ble. The initial conditions for both u (at t — 0) and b (at t = t{) 
correspond to power-law energy spectra with a slope -5/3 cut 
by the equivalent of the dissipation during the duration Sty. 

Each model is run with a wide range of parameters v and r\ 
from 10~ 12 to 1 (6281 independent runs of each model). Such a 
range of parameters is only made possible thanks to the simpli- 
fications operated in the shell models. 
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Figure 1. Growth rate 7 as a function of the Reynolds num- 
bers Re and Rm (note that the color scale is different for positive 
and negative values of y), in the case of the weakly nonlocal 
Sabra model (a - -5/2). The stability curve, corresponding to 
the y = level line, is overplotted. The dashed rectangle cor- 
responds to the range of parameters explored by the direct nu- 
merical simulations of Iskak ov et al.| ([2007); the green triangles 
indicate the initial Reynolds numbers used in Fig.|4j the dotted 
line is where Pm = 1 . 



Every ten time steps, we compute the kinetic energy E u = 
\ Yun \ u n\ 2 , from which we evaluate the hydrodynamic and mag- 
netic Reynolds numbers by Re = 2n y/2E u /kov and Rm = 
2n -yJlEulkoq. We also evaluate the growth rate y of the mag- 
netic energy E/, = I £ l^«l 2 from the local slope of In Et,(t). We 
stop this analysis when £/, becomes greater than lO -4 ^,,, i.e., 
restricting ourselves to the kinematic regime of the dynamo. In 
this way, we get a large set of (Re, Rm, y) triplets, from which 
statistics can be evaluated. 



3. Results 

Growth rate as a function of Reynold numbers. In Fig.[T] we 
plot the average of y in bins of (Re, Rm) for the weakly nonlocal 
model (a = -5/2). We see that the dynamo is effective (7 > 0) 
for higher values of Rm. The level line 7 = represents the 
stability curve Rm f (Re) of the dynamo; a striking feature is that 
it is almost independent of Re and that it seems to have a finite 
limit when Re — > 00. This is a feature that has been noted by 
Iska kov et al.| ( 2007] l up to Re « 3600, and here we confirm this 
result up to values of Re * 10 12 . In addition, at low Re, we also 
recover the decrease in the stability curve. 

The three different models (see also Fig.|5]and|6]online) have 
the same general behavior and the numerical results that we ob- 
tain for Rm c , summarized in Table [T] are almost indistinguish- 
able. One can note however that the stability curves Rm c (Re) are 
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Figure 2. Stability curves Rm c (Re) for the local (plain line), 
weakly nonlocal (dashes), and strongly nonlocal (dots) models, 
and from Iska kov et al.| ([2007 ) (diamonds). 

Table 1. Numerical values obtained for the three models: Re m ; n 
is the smallest Re in Fig. [T] /3 is the slope of the power-law fit 
to /(Rm) = lim Re ^ M y, ft is the slope of the power-law fit to 
/'(Re) = lirriRm^co y, ft is the exponent of the first structure func- 
tion for the velocity field, and y8 p henom. is the expected value for 
ft and ft' deduced from the small-scale dynamo phenomenology. 
Uncertainties in ft, ft' , and ft are 3-cr uncertainties in the fits 
corresponding to these parameters. 



Local GOY 
Local 



Sabra, a = -5/2 
Weak nonloc. 



Sabra, a = - 1 
Strong nonloc. 



lim logRm t . 

Re— >oo 

logRm^Re,™) 
ft 
P 

(ft) 



2.4 ±0.2 



1.5: 

0.424 : 
0.422 : 
0.354: 
0.477 : 



0.3 

0.024 

0.032 

0.067 

0.073 



2.3: 

1.6: 
0.418: 

0.379 : 



0.2 

0.3 

0.032 

0.032 



0.350 ±0.067 
0.481 ±0.073 



2.2 ±0.2 

1.6 ±0.3 
0.379 ± 0.024 
0.336 ± 0.027 
0.383 ± 0.067 
0.446 ± 0.070 




Figure 3. Logarithm of the growth rate y (when positive) as a 
function of the magnetic Prandtl number Pm for different mag- 
netic Reynolds numbers Rm (curves are labeled with log Rm), 
for the weakly nonlocal model. Inset: the limit y^, of y for 
Re — > oo and a power-law fit (plain line). 



shifted upwards (Fig. |2| and that the limit of Rm c as Re — > oo 
slightly increases (Table [T]) when the range of nonlocal interac- 
tions decreases. This implies that nonlocal interactions increase 
the efficiency of the dynamo, even at small Pm where these non- 
local interactions are expected to be weak. 

Growth rates at a given kinetic or magnetic Reynold number. 

The curves y(Pm) for different Rm are plotted in Fig. [3} the plot 
variables are the same as in Fig. 1 of Iskakov et al. (2007). These 
curves correspond to horizontal cuts in Fig.[l] with inverted ab- 
scissa and with a shift by log Rm to the left. They confirm two 
results of Iskakov et al. (2007 ) for the growth rate at a given Rm, 



again over a much wider range of parameters: (1) for a given 
Rm, the growth rate y reaches a limij^] /(Rm) = lim y(Re, Rm) 
for Re — > oo (i.e., Pm — > 0, the saturation of y starting for Pm 
slightly greater than 1); and (2) for any sufficiently large mag- 
netic Reynolds number (Rm > 250), this limit /(Rm) is positive. 

Similarly, vertical cuts in Fig.[T](not shown) demonstrate that 
for a given Reynolds number Re, the growth rate y reaches a 
limit /'(Re) = limy(Re, Rm) for Rm — > oo, i.e., for Pm — > oo. 
Again, the saturation of y starts for Pm slightly greater than 1 ; 
this asymmetry with respect to Pm = 1 can be seen in Fig. [T] as 
a shift between the dotted line (Pm = 1) and the "crest line" of 
the representation of y(Re, Rm), and to our knowledge this shift 
has not been noted before. 



Asymptotic behavior of the limits of the growth rate. We deter- 
mine the limit /(Rm) as defined before by taking the average of 
y on the three leftmost data points of each curve of Fig. [3] with 
the condition that a limit has been reached. The inset of Fig. [3] 
shows that /(Rm) is a power law of Rm: its slope ft is obtained 
by a linear fit log / ~ ft log Rm, and the results for all three mod- 
els are given in Table[T] Similarly, /'(Re) is a power law of Re, 
= and its slope ft' (also given in Table [TJ is obtained by the linear 
fit log/' -yS'logRe. 

Values of ft and ft' are in the range [0.33,0.43]. These val- 
ues for ft and ft', which could not be determined by the direct 
numerical simulations of Isk akov et al.| ( |2007| ), seem to favor a 
dynamo driven by small-scale motions (an exponent 1 /2 is ex- 
pected) over an outer-scale, or mean-field, dynamo (an exponent 
of is expected). Assuming that \u(€)\ ~ (i.e., ft is the expo- 
nent of the first structure function of the velocity field), a small- 
scale dynamo phenomenology (e.g., Stepanov & Plunianf2 008) 
can indeed be summarized as follows: 

- For Pm <k 1, the scale on which the magnetic field grows 
the fastest is the resistive scale t n ~ Rm _1 ^ 1+ ^\ which 
lies within the inertial range of the velocity spectrum. The 
growth rate is the inverse turnover time at this scale, i.e., 
y = u t Jl n ~ Rm (| - fl)/(1+fl) , giving/? = (1 - ft)/(l + ft). 

- Similarly, for Pm » 1, the scale on which the magnetic field 
grows the fastest is the viscous scale t v ~ Re -1 ^ 1+ ^, where 
resistive dissipation is negligible. The growth rate is then 
y = u t Jt v ~ Re<W.)/CH*) t giving ft' = (1 - ft)/(l + ft) = ft. 

With no intermittency, ft = 1/3 and we recover /? = ft' = 1/2 
for the small-scale dynamo, as menti oned before. With hydr ody- 
namic turbulence intermittency^] the | She & Leveque" ( 1994 1 phe- 
nomenology provides the value ft = l/9+2(l-(2/3) 1/J ) « 0.364 
and then the small-scale dynamo phenomenology gives ft = ft' 



0.466, in agreement with the numerical result of Stepanov & 
|Plunian| ( [2dU8l ). 

In our simulations, the values we obtain for ft (shown in 
Table [T] and computed from ten independent runs with v = rj = 
10~ 10 ) are closer to the She & Leveque| 19*94 value than to one- 
third (no intermittency), although both are within 3-cr error bars. 
The small-scale dynamo phenomenology then yields the values 
for j6 and ft' given as j0 p henom.(ft) in Table [T] This prediction 
and the simulation values are lower in the case of the strongly 
nonlocal model, as a consequence of the higher value of ft in 
this case. 



This limit is noted y IX ,(Rm) in 



Iskakov et al. IJ2OO7], but we need to 



use a different notation because of the other limit that we introduce. 
2 A MHD model such as 



Politano & Pouquet 



19951 would not be 



relevant to the kinematic regime of the dynamo (low magnetic field). 
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logi/=-3, logi)=3— 10, At=0.5 



ogl/=-10, logr; = -10, At=0.02 




Figure 4. Evolution of the spectra (compensated by k 5 ^ 3 ) of ki- 
netic energy (black) and magnetic energy (blue), for different 
sets of initial Reynolds numbers (Re, Rm) (shown as green tri- 
angles in Fig. [lj, for the weakly nonlocal Sabra model. The 
spectra are averaged over ten independent runs and are plotted at 
time intervals shown as At in the plot titles. The different shades 
are a guide for understanding the direction of time (from black 
to gray and from dark to light blue, cycling every five spectra). 



The values for /3 and are systematically lower than the 
ySphenom. predictions, but they are still mostly consistent with 
them. They are also lower than the numerical result of Stepanov 



& Plunian] (|2008|l. The difference from the prediction might be 



caused, at least for the nonlocal models, by a contribution of the 
outer-scale dynamo to the growth of the magnetic field, while 
the diff erence from the numerical results of Stepanov & Plunian 
(2008) may come from the different averaging proces^] 



Magnetic energy spectra. The difference, however small, be- 
tween local and nonlocal models is surprising if one considers 
that nonlocal interactions in wavenumber space are expected to 
become important at large Pm. However, for all values of Pm, 
the evolution of magnetic spectra is consistent with mainly lo- 
cal transfers from kinetic to magnetic energy, and so does not 
require important nonlocal transfers: 



L. 



For Pm = 1 (v = t] = 10~ 3 or 10~ 10 in Fig. [4|, magnetic 
energy grows first on scales corresponding to the end of the 
inertial range of the kinetic energy spectrum, close to the dis- 
sipation range, at a speed depending on Reynolds numbers. 
For Pm <s; 1 (v = 10~ 10 and 77 = 10~ 3 in Fig. the 
wavenumber range where magnetic energy can grow is lim- 
ited by the magnetic diffusivity scale. 
For Pm » 1 (v = 10~ 3 and r\ = 10~ 10 in Fig. [4]), the range 
where magnetic energy can grow is limited by the viscous 
scale, showing that transfers from the kinetic to the mag- 
netic energy are still mainly local (despite the inclusion of 
nonlocal terms in the model). 



We note that the scales on which the magnetic energy grows 
most are consistent, at both small and large Pm, with the scales 
used above to compute the magnetic energy growth rate in the 
small-scale dynamo phenomenology. 



3 At each time step, these authors start several simulations from the 
same initial fields, which are the average of the fields computed in at the 
previous time step, while we compute the growth rates from completely 
independent runs of the simulations and compute the average growth 
rate in the end. 



4. Conclusion 

We have computed the growth rate y of the magnetic field in 
the kinematic regime of a dynamo as a function of the Reynolds 
numbers Re and Rm, with, thanks to shell models, a much wider 
range of parameters than previous studies using direct numeri- 
cal simulations. This wider parameter range is of astrophysical 
relevance and brings a new perspective to results from direct nu- 
merical simulations; it allows us to answer some important out- 
standing questions about the kinematic regime of the dynamo, 
assuming of course that results from shell models remain valid 
for the general MHD equations. 

We confirm that the critical magnetic Reynolds number 
Rm c tends to a finite value at large Reynolds numbers Re. 
Furthermore, the growth rate y tends to a finite value /(Rm) when 
Re tends to infinity, and we find a scaling /(Rm) ~ Rrrr 8 with 
f3 « 0.4. Similarly, y tends to a finite value /'(Re) ~ ReP when 
Rm tends to infinity, with {}' » 0.4. Both limits /(Rm) and /'(Re) 
are attained for Rm > Re. 

These scalings can be explained by an intermittent, mainly 
small-scale dynamo. Furthermore, our results imply that nonlo- 
cal interactions (in Fourier space) play a role in the kinematic 
dynamo, although a limited one. This behavior is consistent 
with results on the locality of nonlinear interactions in MHD 
systems: in Alexakis et al.|p005 ), Mi ninni et al.|(|2005) (direct 



numerical simulations), and Plunian & Stepanov (2007} (shell 
models), nonlocal interactions are significant but are mostly con- 
fined to a wavenumber range that is relatively small compared to 
the wavenumber separation factor between shells that interact in 
shell models. 
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Figure 5. (online only) Same as Fig.[T|for the local GOY model. 
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Figure 6. (online only) Same as Fig. [l]for the strongly nonlocal 
(a = -1) Sabra model. 



